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Abstract. PHENIX has measured J/-0 production at backward, central and forward 
rapidities both in p+p and d+Au collisions at ^/snn = 200 GeV during the third run 
of the RHIC collider. From the p+p collisions, we measure the total J/V' cross section. 
From the d+Au collisions, we compute the nuclear modification factor versus transverse 
momentum, rapidity and ccntrality. All results presented here are preliminary. 



1. Introduction 

In nuclear collisions at 200 GeV, J/ips are mainly produced through gluon fusion and 
thus probe the gluon structure function and its modification in nuclei. J/ips are also 
predicted to be very sensitive to the creation of a quark-gluon plasma and are suppressed 
in Pb+Pb collisions at lower energy p. Measuring their yield in d+Au collisions is 
crucial to understand cold nuclear matter effects such as gluon shadowing, transverse 
momentum broadening, parton energy loss or absorption, as well as to set the baseline 
for future Au+Au measurements. 

The PHENIX experiment is able to measure J/ips through their dilepton decay 
in four spectrometers: two central arms covering the mid-rapidity region of |?7| < 0.35 
and twice tt/2 in azimuth and two forward muon arms covering the full azimuth and 
1.2 < \ri\ < 2.4 in pseudorapidity. Electrons are identified in the central arms by their 
Cerenkov rings and by matching the momentum of charged particles reconstructed in 
drift chambers with the energy deposited in an electromagnetic calorimeter. Muons 
are selected by an absorber and identified by the depth they reach in a succession of 
proportional counters staggered with steel walls. The vertex and the event centrality 
are measured by beam-beam counters lying at 3 < |?7| < 3.9. 

In this presentation, we report on J/ip measurements performed during the third 
RHIC run, with 2.74 nb~^ d+Au and 350 nb^^ P+P collisions. For this analysis, we 
use about 300 J/ip from d+Au and 100 from p+p observed in the central arms, and 
1400 and 420 in the muon arms. Previous data samples have already been published 

§ For the full PHENIX author list and acknowledgments, see Appendix " Collaborations" of this volume. 
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Figure 1. J/ip differential 
cross section, multiplied by 
the dilepton branching ratio, 
versus rapidity. 
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Figure 2. Ratio between 
d+Au and p+p J/V" differen- 
tial cross sections, divided by 
2ccccc X f 97, versus rapidity. 



by PHENIX and the analysis presented here is very similar to the one presented in j2] 
(150 nb~^ P+p) and [3] (24 /ib~^ Au+Au). All results are preliminary. 

2. J/ip cross section in proton+proton 

On figure m we show the measured differential cross section versus rapidity, multiplied 
by the dilepton branching ratio. The dielectron measurement is shown as a single 
point at mid-rapidity. The data from each dimuon spectrometer has been split in 
two rapidity bins. Sohd error bars represent statistical and point to point systematic 
uncertainties. The dashed error bars stand for the systematic uncertainties common to 
one spectrometer. An additional 12.3% global error bar is not displayed. 

A fit to a shape generated with PYTHIA using the GRV94HO parton distribution 
is performed and gives a total cross section, multiplied by the dilepton branching ratio 
of 5.9%, equal to: 

BR X a//^ = 159 nb ± 8.5%(fit) ± 12.3%(abs) (1) 

where the first uncertainty comes from the fit and thus includes both the statistical and 
point-to-point systematics. The second uncertainty accounts for absolute systematic 
errors. Variations in the parton distribution functions used to determine the shape was 
demonstrated to be less than 3% |2] and was neglected here. 

3. J/tp production in deuteron+gold 

In d+Au collisions at 200 GeV, J/ips in our three rapidity ranges probe the following 
momentum fractions x of gluons in the gold nucleus (neglecting the emitted gluon): 
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0.05 to 0.14 (backward, negative rapidity, gold-going side), 0.011 to 0.022 (midrapidity) 
and 0.0014 to 0.0047 (forward, positive rapidity, deuteron-going side). 

3.1. Rapidity dependence 

The ratio between the J/ ip yields observed in d+Au collisions and p+p collisions, divided 
by 2 X 197 is shown on figure 121 Solid error bars represent statistical and point to point 
systematic uncertainties. The dashed error bars stand for the systematic uncertainties 
common to one spectrometer. An additional 13.4% global error bar is not displayed. 

While this ratio is close to unity at backward rapidity, it is significantly lower at 
forward rapidity, where parton distribution are expected to be shadowed in a heavy 
nucleus. Theoretical predictions from Vogt |lj and Kopeliovich jH] are displayed on the 
figure for comparison. 

3.2. Transverse momentum dependence 

The J/tp pt distributions {da / dpx) / {^tipt) have been fitted to the traditional (A(l + 
{Pt/BY)~^ function. The average p^ resulting from the fit is 4.47±0.25 and 3.99±0.25 
for d+Au collisions at backward and forward rapidity, respectively, to be compared with 
2.70 ± 0.24 for p+p colhsions. 

Figure El shows the a parameter defined as adA = cTpp x (2A)". Solid error bars 
represent statistical and point to point systematic uncertainties. An additional 2.2% 
global error bar is not displayed. Results from the ^/snn = 38 GeV lower energy 
experiment E866/NuSea |Ej are displayed for comparison. The amounts of pt broadening 
are comparable at both energies. 

3.3. Centrality dependence 

For d+Au collisions, we measure centrality by counting the charge deposited in the 
beam-beam counter facing the gold beam. Centrality is then related to the average 
number of nucleon+nucleon collisions < Ncou > through a Glauber computation. The 
ratio RdA of J/ip yields in d+Au and p+p normalized by the < Ncou > is shown on 
figure m for four centrality classes corresponding to < Ncou >= 3.2 + 0.3, 6.9 + 0.6, 
10.4 + 0.7 and 15.0 + 1.0, and for minimum bias collisions corresponding to an average 
value of Ncoii equal to 8.4 + 0.7. Solid error bars represent statistical and point to point 
systematics and the dashed error bars are common to one spectrometer. An additional 
13.4% global error bar is not displayed. 

At forward rapidity (low x values in the gold nucleus), no strong dependence is 
observed, while a strong enhancement from peripheral to central collisions is observed 
at backward rapidity. The theoretical curves on figurelHstand for corresponding amounts 
of inhomogeneous shadowing and antishadowing at these rapidities from |3]. The 
predictions are qualitatively consistent with the data at forward rapidity, while they 
do not show the observed steep rising shape at backward rapidity. For now, we give no 
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Figure 4. Nuclear modifica- 
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interpretation to this unexpected behaviour, but speculate that it is related to the fact 
that these J ftps are closer to the gold nucleus rest frame. 

4. Conclusions 

During run 3, the PHENIX experiment was able to see nuclear effects affecting J/ip 
production in d-|-Au collisions. The rapidity dependence is consistent with shadowing 
of gluon distribution functions and the px differential cross section shapes exhibit 
broadening. However, given the limited statistics, it is difficult to quantify the level 
of J/ ip absorption or to distinguish between various theoretical models of shadowing. A 
surprising rising shape was seen in the centrality dependence at backward rapidity and 
waits for interpretation. 
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